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Natural  gas  utilities  supply  about  a  quarter  of  the  energy  needs  of  the  United  States.  From  wellhead  to  consumer, 
operations  are  governed  by  an  astounding  diversity  of  purchase,  transport,  and  storage  contract  agreements  which  prepare 
a  complex  physical  distribution  system  to  meet  future  demands  no  more  predictable  than  next  year’s  weather.  We  present 
a  decision  support  system  based  on  a  highly  detailed  optimization  model  used  by  utilities  to  plan  operations  which 
minimize  cost  while  satisfying  regulatory  agencies.  Applications  at  Southwest  Gas  Corporation  are  presented  along  with 
a  case  study  at  Questar  Pipeline  Corporation. 


“But  thou,  contracted  to  thine  own  bright  eyes, 
Feed’st  thy  light’s  flame  with  self-substantial  fuel” 
William  Shakespeare,  First  Sonnet 


Natural  gas  utilities  purchase  gas  from  suppliers, 
arrange  for  its  transport  and  storage,  and 
deliver  it  to  most  end  users  on  an  as-needed  basis.  Gas 
is  distributed  from  the  wellhead  to  consumers  in 
extensive  systems  of  pipelines  and  storage  facilities. 
The  pipelines  can  be  utilized  at  relatively  high  capacity 
even  during  low  demand  periods  because  gas  can  be 
injected  into  underground  storage  facilities  or  lique¬ 
fied  and  stored  in  tanks  for  periods  of  high  demand. 
(Some  gas  is  stored  in  above-ground  gas  holders,  but 
these  are  minor  volumes  used  to  buffer  hour-to-hour 
load  swings.) 

In  1988,  the  United  States  used  about  19  quadrillion 
British  Thermal  Units  (1.9  x  10 16  BTUs1)  of  natural 
gas,  almost  a  quarter  of  the  total  energy  consumption 
in  the  country  (American  Gas  Association  1989b). 


Natural  gas  utilities  supplied  48  million  residential 
customers,  as  well  as  5  million  commerical,  industrial 
and  other  end  users.  The  average  cost  to  end  users 
was  about  $5  per  MMBTU — 10%  less  than  oil,  and 
less  than  a  quarter  the  cost  of  electricity.  Most  of  this 
natural  gas  was  produced  domestically.  Natural  gas 
produces  energy — by  combustion,  or  via  catalysis  and 
indirect  chemical  conversion — with  less  undesirable 
by-products  than  competing  conventional  fuels:  car¬ 
bon  dioxide,  nitrous  oxides,  sulfur  oxides,  and  com¬ 
plex  hydrocarbon  byproducts  are  minimized  (Burnett 
and  Ban  1989). 

Structural  and  regulatory  changes  in  the  natural  gas 
industry  have  increased  the  importance  of  systematic 
analysis  of  supply  options  by  natural  gas  purchasers. 
As  a  result  of  deregulation,  gas  purchasers,  including 
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gas  pipeline  companies,  local  distribution  companies 
(LDCs  or  utilities),  and  large  end  users  are  able  to  take 
advantage  of  new  supply  alternatives.  Developing  an 
optimal  gas  supply  strategy  in  this  new  market  requires 
the  evaluation  of  a  wide  range  of  purchase,  transpor¬ 
tation  and  gas  storage  options  over  a  multiyear  time 
horizon  (Stewart  1987,  American  Gas  Association 
1989a). 

We  describe  a  decision  support  system  called  Gas 
Contract  Analyzer  which  employs  a  linear  program¬ 
ming  model  to  assist  gas  supply  planners  in  making 
both  strategic  and  operational  gas  supply  and  trans¬ 
portation  decisions.  This  system  has  been  used  suc¬ 
cessfully  to  develop  both  short-term  and  long-term 
gas  supply  strategies  by  a  number  of  pipeline  and 
natural  gas  distribution  companies. 

1.  NATURAL  GAS  INDUSTRY  STRUCTURE 

The  natural  gas  industry  consists  of  three  distinct 
segments: 

Producers,  both  major  oil  companies  and  smaller 
independents,  are  primarily  involved  in  exploration 
for,  development  of,  and  production  from  gas  reserves; 
gas  pipeline  companies  aggregate  gas  in  areas  of  pro¬ 
duction  and  transport  the  gas  to  areas  of  consumption; 
LDCs  provide  natural  gas  service  to  end  users  in  a 
particular  geographic  region. 

Despite  recent  moves  toward  deregulation,  govern¬ 
ment  regulatory  bodies  still  have  a  great  deal  to  say 
about  how  natural  gas  is  bought  and  sold.  Although 
production  itself  is  not  regulated,  producer  sales  of 
natural  gas  have  been  subject  to  federal  price  regula¬ 
tion  since  1954.  The  activities  of  interstate  gas  pipeline 
companies  are  regulated  by  the  Federal  Energy 
Regulatory  Commission  (FERC)  under  the  Natural 
Gas  Act  of  1938.  Intrastate  pipelines  and  LDCs  come 
under  the  regulatory  authority  of  state  and  local  public 
utility  commissions. 

LDCs  have  traditionally  purchased  the  majority  of 
their  gas  supplies  from  interstate  pipeline  companies 
under  long-term  contracts.  Contract  terms  can  extend 
twenty  years  or  longer.  These  contracts  allow  the  LDC 
to  purchase  up  to  a  specified  contract  demand  quan¬ 
tity  each  day,  and  up  to  a  maximum  seasonal  quantity 
of  gas  each  year.  Under  these  contracts  the  LDC  pays 
a  monthly  demand  charge  based  on  the  daily  and 
seasonal  contract  quantities,  and  a  commodity  charge 
for  each  unit  of  gas  purchased.  Sales  rates  and  terms 
of  service  are  established  in  pipeline  tariffs  approved 
by  the  FERC. 


An  LDC,  as  a  public  utility  company,  has  a  legal 
obligation  to  provide  natural  gas  service  to  residential, 
commerical,  and  industrial  customers  within  its 
service  territory.  Service  to  larger  commercial  and 
industrial  customers  is  provided  on  either  a  firm  or 
interruptible  basis.  Interruptible  sales  customers  gen¬ 
erally  have  the  capability  of  switching  on  short  notice 
between  natural  gas  and  an  alternate  source  of  fuel, 
such  as  propane  or  fuel  oil.  In  exchange  for  a  lower 
delivered  price,  interruptible  customers  face  the  risk 
that  gas  supplies  will  be  curtailed  during  periods  of 
peak  demand. 

LDCs  meet  peak  demand  requirements  through  a 
combination  of  contract  purchases  and  withdrawals 
from  storage.  Gas  is  stored  in  underground  or  liquefied 
natural  gas  facilities  operated  by  the  LDC,  or  through 
storage  service  contracts  with  pipeline  suppliers.  In 
either  case,  the  LDC  incurs  a  cost  related  to  the 
amount  of  storage  capacity  available,  a  cost  related  to 
the  maximum  quantity  that  can  be  withdrawn  from 
storage  on  a  given  day,  and  a  cost  for  each  unit  of  gas 
which  is  injected  into  or  withdrawn  from  the  storage 
facility. 

2.  OPEN  ACCESS 

The  FERC  has  implemented  a  series  of  regulatory 
changes  (beginning  with  Order  436  in  1 985)  to  encour¬ 
age  pipeline  companies  to  unbundle  their  traditional 
merchant  services  and  allow  buyers  to  transport  their 
own  gas  through  the  pipeline  system.  This  open  access 
to  gas  transportation  services  gives  gas  purchasers  the 
opportunity  to  buy  gas  directly  from  nonpipeline 
sources,  including  natural  gas  producers  and  inde¬ 
pendent  marketing  companies,  and  has  contributed 
to  the  development  of  an  active  spot  market  for 
natural  gas. 

As  part  of  this  market  restructuring,  pipeline  cus¬ 
tomers  have  been  given  the  opportunity  to  convert  all 
or  part  of  their  existing  gas  purchase  entitlements  to 
firm  transportation  service.  A  firm  transportation  con¬ 
tract  allows  the  shipper  to  reserve  a  portion  of  the 
pipeline’s  total  delivery  capacity  for  his  own  use.  The 
shipper  pays  a  monthly  demand  charge  based  on  the 
maximum  daily  delivery  quantity  contracted  for,  and 
a  transportation  charge  for  each  unit  of  gas  delivered. 
Additional  pipeline  transportation  service  is  available 
on  an  interruptible  basis.  For  interruptible  transpor¬ 
tation  services  the  shipper  generally  pays  only  for  the 
gas  transported. 

Open  access  has  had  perverse  consequences  for 
some  pipeline  companies  and  LDCs.  Major  industrial 
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customers  of  an  LDC  have  traditionally  contracted 
for  inexpensive  interruptible  gas  from  the  LDC  which 
obtains  the  gas  from  some  pipeline  company.  In  the 
past,  interrupting  these  customers  has  provided  the 
LDC  with  a  handy  method  of  leveling  demand  peaks. 
Under  open  access,  however,  many  of  these  customers 
have  contracted  directly  with  a  producer  to  supply 
gas,  with  a  pipeline  company  being  obligated  to  pro¬ 
vide  interruptible  transportation.  Typically,  the  LDC 
will  have  concurrent  sales  agreements  with  the  indus¬ 
trial  customers  which  stipulate  that  if  the  producer 
has  trouble  meeting  its  commitments,  the  LDC  will 
make  up  the  deficit  at  “best  efforts.”  Thus,  the  indus¬ 
trial  customers  have  become  a  part  of  the  peak-load 
problem  rather  than  contributing  to  its  solution  and, 
in  fact,  in  the  last  few  years  the  proportion  of  inter¬ 
ruptible  demand  as  a  fraction  of  total  demand  has 
decreased  by  one  third. 

Finally,  because  gas  buyers  are  now  able  to  swing 
from  one  supplier  to  another  based  on  price,  many 
gas  contracts  include  pricing  provisions  to  encourage 
buyers  to  maintain  high  purchase  load  factors.  At  the 
pipeline  level,  many  companies  have  proposed  new 
gas  inventory  charges.  These  charges  are  structured 
either  as  a  demand  charge  based  on  the  buyer’s  gas 
purchase  entitlement,  or  as  a  penalty  charge  for  each 
unit  that  total  monthly  or  seasonal  purchases  fall  short 
of  a  contractual  target  quantity.  Direct  sales  contracts 
arranged  with  producers  or  independent  marketers 
often  contain  similar  provisions.  In  the  extreme,  the 
buyer  is  required  to  pay  for  a  specified  quantity  of  gas 
whether  it  is  taken  or  not. 

3.  GAS  SUPPLY  OBJECTIVES 

Purchased  gas  costs  can  account  for  60%  or  more  of 
the  total  costs  of  an  LDC.  The  objective  of  gas  supply 
planning  is  to  minimize  gas  supply  costs  while  main¬ 
taining  sufficient  supply  to  meet  potential  peak 
requirements  and  provide  for  future  growth  in 
demand. 

Over  the  short  term,  this  means  optimally  dispatch¬ 
ing  the  available  gas  supply  to  meet  variable  demand. 
A  gas  dispatch  plan  reconciles  gas  purchases,  storage 
inventory  targets  for  key  points  in  the  injection/ 
withdrawal  cycle,  and  the  expected  level  of  firm  and 
interruptible  sales. 

Over  the  long  term,  the  objective  of  supply  planning 
is  to  construct  an  optimal  portfolio  of  gas  sources 
including  gas  purchases  from  pipelines,  storage,  and 
transportation  of  gas  purchased  directly  from  the  pro¬ 
ducer.  As  the  menu  of  gas  supply  options  has 
increased,  the  problem  of  sorting  through  these 


options  to  develop  a  coherent  gas  supply  strategy  has 
become  much  more  complex.  Among  the  strategic 
issues  that  supply  planners  must  now  address  are  the 
following: 

Contract  Restructuring 

•  How  much  of  the  company’s  existing  pipeline 
purchase  entitlement  should  be  converted  to  firm 
transportation? 

•  Should  the  overall  entitlement  level  be  reduced  or 
increased? 

Direct  Purchases 

•  How  much  gas  should  be  purchased  from  nonpipe¬ 
line  suppliers  under  long-term  contracts? 

•  What  role  should  short-term  spot  purchases  play  in 
the  company’s  supply  portfolio? 

Storage 

•  Should  the  company  develop  or  lease  additional 
storage  capacity? 

•  What  storage  capacity  and  deliverability  levels  are 
best? 

Marketing 

•  What  types  of  gas  suppliers  should  support  each 
customer  market  segment? 

•  How  much  gas  should  the  company  make  available 
for  interruptible  sales? 

4.  PROBLEM  CONSTRUCTS 

We  seek  a  flexible  optimization  system  designed  to 
address  both  contract  restructuring  decisions  and  the 
optimal  utilization  of  gas  supply  sources.  The  key 
elements  of  the  optimization  problem  are  outlined 
below. 

Time  Periods 

The  dispatch  periods  can  be  daily,  weekly,  monthly, 
or  any  aggregation  of  these.  The  planning  horizon  can 
range  from  one  year  to  over  a  decade. 

Gas  Contracts 

All  sources  of  gas  are  modeled  as  contracts.  The  model 
sets  dispatch  flow  for  all  contracts,  each  of  which  is 
assumed  to  have  its  own  commodity  rate,  i.e.,  a  cost 
proportional  to  the  amount  of  gas  purchased.  In  addi¬ 
tion,  the  model  determines  optimal  maximum  daily 
and  maximum  seasonal  purchase  levels  within  pre¬ 
defined  upper  and  lower  bounds.  For  potential  con¬ 
tracts  (as  opposed  to  existing  contracts),  a  D1  charge 
is  incurred  proportional  to  the  maximum  daily 
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purchase  level  and  a  D2  charge  is  similarly  incurred 
for  the  maximum  seasonal  purchase.  Both  existing 
and  potential  contracts  can  also  be  assigned  a  mini¬ 
mum  seasonal  purchase  quantity,  and  a  penalty  rate, 
called  a  deficiency-based  gas  inventory  charge,  to  be 
applied  to  any  purchase  shortfall.  Groups  of  contracts 
may  share  joint  daily  and  seasonal  purchase  limits. 

Storage 

Existing  gas  storage  facilities  or  potential  storage 
service  contracts  have  a  maximum  daily  injection 
volume,  a  maximum  daily  withdrawal  volume,  and  a 
maximum  storage  capacity.  For  potential  contracts, 
within  preset  lower  and  upper  bounds,  the  maximum 
daily  withdrawal  level  can  be  determined  taking  into 
account  the  storage  service  demand  charge  and  the 
maximum  storage  capacity  can  be  determined  taking 
into  account  the  storage  service  capacity  charge. 
Optional  storage  relationships  can  be  used  to  control 
the  rate  of  withdrawals  and/or  injections  as  a  function 
of  gas  inventory.  These  relationships  reflect  the  fact 
that  the  higher  inventory  is  the  easier  it  is  to  withdraw, 
but  the  harder  it  is  to  inject. 

Gas  Transportation 

The  dispatch  of  gas  from  sources  to  dispositions  can 
be  characterized  by  defining  a  network  of  nodes  and 
arcs.  Each  arc  can  be  given  minimum  and  maximum 
daily  flow  values,  a  charge  per  unit  of  throughput 
called  the  transportation  rate,  and  a  loss  percentage 
on  each  unit  transported.  For  potential  transportation 
contracts,  the  maximum  throughput  capacity  is  set  for 
both  firm  and  interruptible  transportation  service 
based  on  costs  called  the  firm  demand  rate  and  the 
interruptible  demand  rate,  respectively. 

Side  Constraints 

Additional  dispatch  constraints  can  be  defined  linking 
flows  across  network  arcs  and  across  time  periods. 
This  feature  is  used  to  model  contractual  limitations 
on  the  sources  used  to  inject  into  storage,  and  to 
segregate  transportation  for  end  users  from  system 
supply  volumes  in  a  commingled  gas  flow. 

Peak-Day  Constructs 

The  model  determines  optimal  gas  flows  for  each 
dispatch  period.  In  addition,  the  sizing  of  supply, 
storage,  and  transportation  must  satisfy  a  parallel  sys¬ 
tem  of  peak-day  constraints.  The  peak-day  constraints 
require  that  total  system  deliverability,  adjusted  for 
the  expected  peak-day  reliability  of  existing  and  poten¬ 
tial  gas  sources,  must  be  sufficient  to  meet  peak-day 


demands.  Separate  peak-day  constraints  can  be  spec¬ 
ified  for  any  aggregate  time  period,  such  as  a  month. 

The  parallel  peak-day  submodel  provides  a  mana- 
gerially  appealing  mechanism  to  account  for  future 
demands  which  cannot  be  forecasted.  Reliability  fac¬ 
tors  express  in  common  parlance  the  probability  that 
agreed  limits  will  actually  be  met  under  peak  condi¬ 
tions:  On  a  peak  day  some  contractual  partners  are 
more  trustworthy  than  others,  and  some  contracts  are 
easier  to  fulfill  than  others.  Using  the  reliability  factor, 
such  quantities  as  maximum  daily  purchase  for  a 
pipeline  contract  are  replaced  by  their  expected  value 
on  the  peak  day.  More  elegant  stochastic  program¬ 
ming  methods  might  require  more  data  than  can  be 
certified  comfortably  and  typically  yield  optimization 
problems  which  are  much  larger  and  much  harder  to 
solve  than  our  corresponding  deterministic  problems. 

Objective  Function 

The  present  value  of  all  costs  is  minimized.  This 
includes  variable  costs  on  all  gas  sources,  storages,  and 
transportation  arcs;  committed  charges  on  all  con¬ 
tracts,  storages,  and  arcs  for  which  a  sizing  decision  is 
made;  and  penalty  costs  associated  with  failure  to  meet 
gas  deliverability  requirements  and  minimum  pur¬ 
chase  levels.  A  slight  overcounting  occurs  in  the  objec¬ 
tive  function  because  peak-day  variables  also  have 
costs  associated  with  them.  However,  such  costs  are 
necessary  to  yield  sensible  peak-day  flow  patterns; 
leaving  peak-day  costs  at  zero  could  lead  to  peak-day 
flows  from  the  most  distant,  most  expensive  supplier 
with  the  greatest  amount  of  transmission  loss. 

We  assume,  for  the  sake  of  tractability,  that  all  costs 
are  linear;  we  ignore  the  fact  that  nonlinear  costs  can 
arise  when  constructing  storage  facilities  and  when 
purchasing  gas.  In  the  first  case,  the  modeler  will 
normally  be  considering  discrete  options,  such  as 
opening  two  caverns  of  a  salt  dome,  or  four,  or  maybe 
six.  Typically,  the  different  options  will  be  modeled 
as  existing  storages  in  separate  scenarios  rather  than 
modeling  the  options  in  a  single  scenario  as  a  storage 
service  contract.  Thus,  construction  costs  are  com¬ 
puted  completely  outside  of  the  model.  The  second 
case  occurs  when  there  are  quantity  discounts  for  gas 
purchase  contracts.  In  this  case,  the  modeler  would 
normally  apply  the  lowest  price  and  determine  if  the 
model  purchases  sufficient  gas  to  achieve  that  dis¬ 
counted  price.  If  it  does,  the  solution  is  optimal.  If 
not,  the  next  higher  price  could  be  used  and  the  model 
rerun  to  see  if  the  appropriate  level  of  purchases  is 
made,  etc.  While  optimality  cannot  be  guaranteed, 
this  technique  seems  to  work  well  in  practice  and 
we  have  not  felt  compelled  to  introduce  the  binary 
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variables  and  additional  constraints  necessary  to 
model  these  quantity  discounts  exactly. 


5.  MODEL  CONSTRUCTS 

Our  model  of  gas  transmission  by  pipeline  and  distri¬ 
bution  by  LDCs  is  deceptively  simple  to  visualize. 
Superficially,  a  single  commodity,  base-demand  nat¬ 
ural  gas,  is  purchased  at  supply  points,  held  in  storage 
facilities,  and  transported  through  a  pipeline  network 
to  meet  the  demands  of  different  customer  classes 
over  many  time  periods.  Losses  resulting  from  con¬ 
sumption  of  gas  in  compressors  and  other  line  losses 
are  modeled  both  on  transportation  links  and  in  stor¬ 
age  facilities. 

However,  this  simple  single-commodity  generalized 
network  flow  model  of  base  demand  fails  to  capture  a 
host  of  essential  details.  For  instance,  providing  for 
peak-day  capacity  requires  a  coupled,  parallel  model 
for  meeting  extreme  demands,  incorporating  “peak- 
day  gas”  flowing  through  a  “peak-day  network.”  In 
addition,  contract  stipulations  may  require  that  the 
gas  provided  be  treated  distinctly  from  other  gas. 
Variables  and  constraints  related  to  setting  the  size 
and  other  parameters  of  purchase,  storage  and  trans¬ 
mission  contracts  complicate  the  model  further.  Thus, 
the  resulting  model  has  multiple  commodities,  com¬ 
plicating  constraints  and  variables,  and  bears  a  distant 
resemblance,  at  best,  to  a  simple  generalized  network 
flow  model. 

Incorporation  of  all  these  complicating  features  into 
a  single,  unified  optimization  model  (linear  program) 
is  a  distinguishing  feature  of  the  work  reported  here. 
The  objective  is  to  minimize  cost,  while  meeting  all 
firm  demands  and  interruptible  demands  with  best 
efforts.  We  assume  that  all  data  are  deterministic 
although  we  discuss  how  uncertainty  in  demand, 
availability  and  price  is  incorporated  in  Section  6. 

The  model  consists  of  a  template  of  numerous 
standard  constructs  which  may  be  selected,  specified, 
and  connected  in  various  ways  to  represent  a  gas 
supply  system.  Below,  we  describe  these  building 
blocks,  first  those  associated  with  the  period-by-period 
aspects  of  the  model,  and  then  those  associated  with 
the  peak-day  aspects.  The  parenthetic  mnemonics  key 
and  a  complete  mathematical  formulation  are  given 
in  the  Appendix. 

Period-By-Period  Constructs 

Spot  Purchase  Node.  A  spot  purchase  node  is  a  sim¬ 
ple  supply  node  with  upper  and  lower  limits  on  the 
amount  of  gas  which  can  be  purchased  in  any  period 


as  specified  by  maximum  and  minimum  daily  pur¬ 
chase  levels  ( MXMP ).  A  cost  proportional  to  the 
amount  purchased  is  incurred. 

Existing  Pipeline  Contract  Node.  In  addition  to  the 
features  of  a  spot  purchase  node,  over  a  specified 
season  {MXMP),  an  upper  bound  on  total  purchases 
for  a  season,  i.e.,  a  maximum  seasonal  purchase,  is 
specified  together  with  a  fraction  of  the  maximum 
seasonal  purchase  which  must  be  purchased  ( MXSP ). 
This  minimum  purchase  may  be  violated  at  a  specified 
cost  per  unit  of  violation. 

Potential  Pipeline  Contract  Node.  A  pipeline  contract 
may  be  fixed  for  the  next  few  years,  and  then  be 
subject  to  renegotiation.  During  the  fixed  period,  its 
specification  is  identical  to  that  of  an  existing  pipeline 
contract  node.  After  that  time,  the  following  char¬ 
acteristics  associated  with  a  potential  contract  are 
determined: 

•  A  maximum  daily  purchase  level  is  chosen  and 
enforced  (MXMM)  and  a  cost  proportional  to  the 
level  chosen  is  incurred. 

•  A  fixed  minimum  daily  purchase  may  be  specified 
(. MNMP ). 

•  A  maximum  seasonal  purchase  level  is  chosen  and 
enforced  (SEAS)  and  a  cost  proportional  to  the  level 
chosen  is  incurred. 

•  Also,  a  minimum  seasonal  purchase  may  be  speci¬ 
fied  as  a  given  fraction  of  the  maximum  seasonal 
purchase  level  being  chosen  ( MNPF ).  This  con¬ 
straint  may  be  violated  at  a  specified  cost  per  unit 
of  violation. 

•  The  maximum  daily  purchase  and  the  maximum 
seasonal  purchase  may  be  related  to  each  other 
(DSLO  and  DSUP). 

Parent-Children  Pipeline  Contract  Nodes.  In  addition 
to  the  above  single-node  contracts,  it  is  possible  to 
specify  contracts  which  cover  several  nodes  in  a 
parent-to-children  relationship.  In  addition  to  the  con¬ 
straints  associated  with  the  individual  nodes,  the  fol¬ 
lowing  constraints  may  be  specified: 

•  A  joint  maximum  daily  purchase  level  (over  the 
parent  and  its  children)  is  chosen  and  enforced 
{JMXM),  while  a  cost  proportional  to  that  level  is 
incurred. 

•  A  joint  maximum  seasonal  purchase  level  is  chosen 
and  enforced  ( JSEA ),  while  a  cost  proportional  to 
that  level  is  incurred. 

•  A  maximum  daily  purchase  level  over  the  children 
in  the  contract  is  enforced  as  a  fraction  of  the  joint 
maximum  daily  purchase  ( JCPM ). 
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•  The  joint  maximum  daily  purchase  and  the  joint 

maximum  seasonal  purchase  may  be  related  to  each 

other  (JDSL  and  JDSU). 

Existing  Storage  Node.  This  is  structurally  modeled 
as  an  injection  node,  an  inventory  node  and  a  with¬ 
drawal  node  with  bypasses  allowed,  i.e.,  gas  may  be 
transferred  through  the  storage  location  without 
actually  entering  and  leaving  storage.  Loss  fractions 
for  both  injections  and  withdrawals  may  be  specified 
along  with  minimum  and  maximum  levels  for  these 
quantities. 

Withdrawals  in  any  period  may  be  limited  to  a 
fraction  of  the  total  available  inventory  plus  an  offset. 
This  approximates  the  phenomenon  that  withdrawals 
from  inventory  become  more  difficult  as  inventory 
levels  drop  ( WINV). 

Injections  in  any  period  may  be  limited  by  an 
inverse  function  of  the  inventory.  This  approximates 
the  phenomenon  that  injections  into  inventory 
become  more  difficult  as  inventory  levels  rise  ( IJIV ). 

Flow  into  the  “injection  node”  equals  the  amount 
of  gas  injected  into  storage  plus  the  amount  bypassed 
( INPT ).  A  cost  proportional  to  the  amount  injected 
is  incurred. 

Inventory  during  any  time  period  equals  total  injec¬ 
tions  less  total  withdrawals  up  to,  but  not  including, 
that  period,  and  is  limited  by  total  inventory  capacity. 
A  minimum  inventory  target  level  may  be  specified 
with  a  per-unit  penalty  for  violating  the  target 
(STOR). 

The  amount  of  gas  leaving  a  storage  facility  equals 
withdrawals  plus  gas  which  bypasses  actual  storage 
(OUTP).  A  cost  proportional  to  the  amount  with¬ 
drawn  is  incurred. 

Potential  Storage  Node.  A  potential  storage  node  has 
the  same  physical  structure  as  an  existing  storage  node 
and  the  same  constraints  dealing  with  maximum  daily 
withdrawals  (WINV),  maximum  daily  injections 
(IJIV),  input  (INPT),  and  output  (OUTP). 

In  contrast,  the  maximum  inventory  level  is  not 
fixed  but  is  chosen  by  the  model  within  specified  limits 
(STOR).  A  cost  proportional  to  this  level  is  incurred. 

A  separate  constraint  enforces  inventory  targets  if 
desired  (INVT). 

The  maximum  daily  withdrawal  is  also  a  variable 
which  must  be  selected  by  the  model  within  specified 
limits  (VMDW).  A  cost  proportional  to  the  maximum 
daily  withdrawal  is  incurred. 

Junction  Node.  A  junction  node  requires  that  inflow 
equal  outflow  (JUNC). 


Demand  Node.  A  demand  node  simply  requires  that 
inflow  less  outflow  must  equal  demand,  except 
that  demand  may  be  shortened  at  a  specified  cost  per 
unit  (DMDM). 

Existing  Transportation  Arc.  This  is  simply  a  link 
connecting  two  nodes  in  the  network  with  minimum 
and  maximum  daily  throughput  constraints.  A  loss 
factor  may  also  be  specified  along  with  a  cost  per  unit 
of  gas  transported. 

Potential  Transportation  Arc.  Two  parallel  links  are 
used  to  model  this  situation,  representing  firm  and 
interruptible  transportation.  These  links  have  the 
properties  of  existing  transportation  arcs  along  with 
the  following  additional  properties: 

•  The  total  flow  over  the  two  links  is  limited  by  a 
joint  capacity  constraint  (JCAP). 

•  A  maximum  daily  throughput  is  chosen  and 
enforced  individually  for  both  firm  and  interruptible 
transportation  (MXMF and  MXMI).  Separate  costs, 
proportional  to  the  maximum  throughputs  are 
incurred. 

•  Interruptible  transportation  may  be  reduced  by  a 
multiplier  ranging  from  0  to  1  in  periods  when  it  is 
partially  or  completely  unavailable. 

Side  Constraints.  The  flow  on  any  number  of  trans¬ 
portation  arcs  may  be  linked  together  over  a  specified 
set  of  time  periods  (ALNK). 

Peak-Day  Constructs 

The  parallel  peak-day  submodel  ensures  that  demands 
can  be  met  on  a  peak  day  as  well  as  on  an  average 
day.  This  submodel  dispatches  gas  on  the  peak  day 
not  in  addition  to  the  period-by-period  requirements 
but  subject  to  the  contract  constraints  described 
above.  Additionally,  the  withdrawals  made  from  stor¬ 
age  on  the  peak  day  are  limited  by  the  inventory  levels 
obtained  at  the  end  of  the  period.  Peak-day  constructs 
are  automatically  generated  when  their  parallel  non¬ 
peak  constructs  are  selected,  although  the  peak-day 
constructs  may  be  limited  to  specific  periods  when 
peaking  capabilities  are  important,  i.e.,  during  the 
winter.  The  individual  constructs  are  described  below. 

Peak-Day  Spot  Purchase  Node.  Total  purchases  are 
limited  by  minimum  and  maximum  daily  purchase 
levels  possibly  adjusted  downward  by  a  “reliability 
factor”  (PMDP).  A  cost  proportional  to  the  amount 
purchased  is  incurred. 

Peak-Day  Existing  Pipeline  Contract  Node.  This  is 
the  same  as  a  peak-day  spot  purchase  node  (PMDP). 
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Peak-Day  Potential  Pipeline  Contract  Node.  Peak 
purchases  are  limited  by  the  maximum  daily  purchase 
level  chosen,  possibly  adjusted  downward  by  a  relia¬ 
bility  factor  ( PMMD ).  A  cost  proportional  to  the 
amount  purchased  is  incurred. 

The  peak  purchases  must  exceed  the  minimum 
daily  purchase  requirement  ( PMND ). 

For  a  joint  (parent-children)  contract,  the  total  pur¬ 
chases  may  not  exceed  the  joint  maximum  daily  pur¬ 
chase  level  chosen,  possibly  adjusted  downward  by  a 
reliability  factor  ( PJMM ). 

For  a  joint  contract,  the  total  purchases  from  the 
children  may  not  exceed  a  given  fraction  of  the  joint 
maximum  daily  purchase  level  chosen  ( PJCC ). 

Peak-Day  Existing  Storage  Node.  Unlike  the  non¬ 
peak  portion  of  the  model,  a  peak-day  existing  storage 
node  is  structurally  a  single  node. 

Outflow  from  the  node  equals  inflow  plus  with¬ 
drawals  ( PSTR ). 

Withdrawals  are  limited  to  a  fraction  of  the  period¬ 
ending  inventory  plus  an  offset  ( PINV ).  A  cost  pro¬ 
portional  to  the  withdrawal  is  incurred. 

In  no  case  may  withdrawals  exceed  inventory 
(. PTIV )  or  the  maximum  daily  withdrawal  which  is 
possibly  adjusted  downward  by  a  reliability  factor. 

Peak-Day  Potential  Storage  Node.  A  potential  stor¬ 
age  node  has  the  constraints  {PSTR),  {PINV)  and 
{PTIV)  as  in  an  existing  storage  node.  However,  the 
fixed  upper  bound  on  withdrawals  is  replaced  by  a 
variable  maximum  daily  withdrawal  possibly  adjusted 
downward  by  a  reliability  factor  {PVMW). 

Peak-Day  Junction  Node.  Flows  into  the  node  equal 
flows  out  of  the  node  on  the  peak  day  {PJNC). 

Peak-Day  Demand  Node.  Demand  on  the  peak  day 
equals  gas  flowing  into  the  node  less  gas  flowing  out 
of  the  node.  Demand  can  be  shorted  at  a  specified 
cost  per  unit  {PDMD). 

Peak-Day  Existing  Transportation  Arc.  This  is  sim¬ 
ply  a  single  link  connecting  two  nodes  in  the  peak- 
day  submodel.  Row  is  limited  by  the  minimum  daily 
throughput  and  the  maximum  daily  throughput,  pos¬ 
sibly  adjusted  downward  by  a  reliability  factor.  A  cost 
per  unit  of  gas  transported  is  incurred. 

Peak-Day  Potential  Transportation  Arc.  This  is  rep¬ 
resented  as  two  links,  in  parallel,  connecting  two  nodes 
in  the  peak-day  portion  of  the  model.  Costs  propor¬ 
tional  to  the  amount  of  flow  on  each  link  are  incurred. 


The  sum  of  the  flow  on  the  links  is  limited  by  their 
joint  capacity  {PTCP). 

The  peak-day  flow  on  the  firm  transportation  link 
is  limited  by  the  maximum  daily  (firm)  throughput 
chosen  by  the  model.  This  may  be  adjusted  downward 
by  a  reliability  factor  {PMXF). 

The  peak-day  flow  on  the  interruptible  transporta¬ 
tion  link  is  limited  by  the  maximum  daily  (interrup¬ 
tible)  throughput  chosen  by  the  model.  This  may  be 
adjusted  downward  by  a  reliability  factor  which  is 
distinct  from  the  “firm”  reliability  factor  {PMXI). 

6.  DEALING  WITH  UNCERTAINTY 

While  Contract  Analyzer  handles  variability  of 
demand  and  peak-day  surge  it  is  essentially  a  deter¬ 
ministic  model.  It  is  meant  to  be  used  for  analyzing  a 
number  of  future  demand/price  scenarios  posited  by 
an  LDC.  Typically,  a  number  of  scenarios  are  devel¬ 
oped  and  run  under  Contract  Analyzer  and  solutions 
are  compared.  If  there  is  little  difference  among  the 
scenarios  in,  say,  the  purchase  level  of  a  particular 
contract,  the  analyst  can  be  fairly  certain  of  the  correct 
purchase  level.  If  instead,  significant  variability  arises, 
the  analyst  can  determine  a  composite  purchasing 
strategy,  fix  it  and  then  run  the  various  scenarios 
against  that  strategy  to  see  how  solutions  vary.  Sensi¬ 
tivity  analyses  on  the  strategy  can  then  be  carried  out 
if  desired.  (See  Wagner  1969,  Chapter  16  for  a  classic 
discussion  of  this  type  of  analysis.) 

While  it  would  be  possible  to  create  a  two-  or 
multistage  stochastic  programming  model  with 
recourse  (e.g.,  Walkup  and  Wets  1967)  which  would 
simultaneously  incorporate  different  scenarios,  the 
size  of  such  a  problem  (it  would  grow  at  least  linearly 
with  the  number  of  scenarios  included)  would  make 
its  solution  very  difficult.  Scenario  analysis  gives  good 
results  in  practice  and  its  use  is  well  accepted  by  LDCs 
and  by  their  regulators. 

Note  that  a  stochastic  model,  called  the  “Contract 
Mix  Model”  (Fancher,  Wilson  and  Mueller  1985)  can 
be  used  for  analysis  of  contract  purchases  under  uncer¬ 
tain  demand  and  prices.  While  this  model  may  be 
useful  in  some  circumstances  it  has  a  number  of  severe 
limitations.  In  particular,  the  model  does  not  consider 
storage  and  transportation  contracts,  it  cannot  model 
a  gas  distribution  network,  it  works  only  at  a  yearly 
level  of  detail,  and  requires  strong  assumptions  about 
the  independence  of  demands  and  prices. 

To  develop  scenarios  for  analysis,  LDCs  must  try 
to  predict  likely  levels  of  future  demand  and  costs. 
Predicting  demand  for  natural  gas  is  neither  easier  nor 


more  reliable  than  long-range  climatological  and 
econometric  forecasts:  Demand  is  dominated  by 
domestic  space  heating  and  industrial  applications. 
However,  utilities  normally  have  available  climatolog¬ 
ical  data  for  a  past  “normal”  year,  a  cold  year  and  a 
very  cold  year.  Using  these  data,  along  with  historical 
demand  for  those  years  and  the  best  econometric 
models  available,  an  LDC  can  develop  reasonable 
scenarios  for  future  demands. 

Price  scenarios  are  more  difficult  to  develop  in  a 
logical  fashion.  Natural  gas  prices  in  the  future  will 
depend  on  demand,  gas  availability,  availability  and 
price  of  alternate  energy  sources,  and  politics.  LDCs 
have  used  statistical  forecasting  models  for  predicting 
future  prices,  but  these  are  usually  tempered  by  man¬ 
agerial  judgment  to  consider  such  things  as  spikes  in 
prices  caused  by  unforecastable  economic  and  politi¬ 
cal  events.  At  least  one  LDC  is  planning  on  developing 
price  scenarios  using  a  Delphi  method  to  incorporate 
the  judgment  of  its  best  analysts  and  managers. 

7.  DECISION  SUPPORT  SYSTEM  SHELL 

Our  system  provides  a  user-friendly  interface  for  the 
optimization  model.  Data  entry,  scenario  manage¬ 
ment,  report  writing,  and  ad  hoc  analysis  of  results 
are  supported  via  full-screen  interfaces  and  command 
menus. 

Data  entry  is  supported  by  full-screen  templates 
which  prompt  the  required  data  and  automatically 
perform  range  and  consistency  checks.  This  reduces 
false  starts  when  data  are  being  debugged.  Data  are 
organized  in  a  compact  form  so  that  changes  can  be 
made  without  extensive  regeneration. 

Scenarios  are  generated  by  pointing  to,  rather  than 
copying  data.  Thus,  alternate  excursions  are  easy  to 
assemble  via  command  macros  to  set  conditions  for 
each  run.  By  contrast,  traditional  data  management 
schemes  generate  multiple  copies  of  input  data  sets, 
one  for  each  scenario.  Each  of  the  data  sets  must  be 
individually  edited,  validated,  and  debugged  making 
coordination  of  multiple  scenarios  difficult.  For 
instance,  if  one  fundamental  constant  is  changed, 
traditional  systems  require  the  editing  of  the  data  for 
each  of  the  multiple  scenarios.  An  additional  advan¬ 
tage  of  our  system  is  that  the  user  is  not  expected 
to  master  editors,  job  control  language  or  other 
system-specific  details.  In  fact,  the  user  sees  very  few 
differences  between  implementations  on  different 
computers  or  operating  systems. 

Report  writing  is  easy.  Standard  reports  display 
most  of  the  interesting  details.  Report  modifications 
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are  made  with  simple  English  commands  to  format 
output  displays.  These  commands  are  interpreted  and 
then  applied  to  a  standard  output  data  template  to 
render  the  required  report. 

8.  APPLICATIONS  AT  SOUTHWEST  GAS 
CORPORATION 

To  date,  Gas  Contract  Analyzer  has  been  installed  at 
eight  LDCs  around  the  country.  The  use  of  this  model 
at  Southwest  Gas  Corporation  is  typical. 

Southwest  Gas  Corporation  operates  gas  transmis¬ 
sion  and  distribution  facilities  in  the  states  of  Arizona, 
Nevada,  and  California.  In  1 989,  Southwest  Gas  deliv¬ 
ered  133,000  MMCF  to  797,000  customers.  Sales 
comprised  68%  of  this  volume,  while  transportation 
of  customer-owned  gas  comprised  the  remaining  32%. 
Revenues  from  gas  operations  totaled  $570  million. 

Southwest’s  gas  distribution  system,  depicted  in 
Figure  1,  is  divided  into  the  Northern  System, 
covering  northern  Nevada  and  a  small  portion  of 
California,  and  the  Southern  System,  which  serves 
southern  Nevada,  Arizona,  and  extreme  southeastern 
California.  Each  system  is  connected  to  a  single  inter¬ 
state  pipeline.  Northwest  Pipeline  System  serves  the 
Northern  distribution  system  through  Southwest’s 
interstate  gas  transmission  subsidiary,  Paiute  Pipeline 
Company.  El  Paso  Natural  Gas  (EPNG)  serves  the 
Southern  System.  Peak  deliverability  on  the  Northern 
System  is  supplemented  by  a  liquefied  natural  gas 
storage  facility  and  a  liquid  petroleum  gas  plant.  There 
is  no  connection  between  the  two  Southwest  Gas 
distribution  systems  at  the  present  time. 

With  open  access  transportation,  Southwest  Gas  has 
been  able  to  purchase  an  increasing  portion  of  its  gas 
from  nonpipeline  sources.  Early  in  1989,  Southwest 
Gas  converted  100%  of  its  sales  entitlement  on 
Northwest  Pipeline  to  firm  transportation.  Spot  mar¬ 
ket  purchases  now  account  for  approximately  60%  of 
total  gas  purchases  system-wide. 

Southwest  Gas  is  currently  restructuring  its  gas  sup¬ 
ply  arrangements  for  the  Southern  System.  The  gas 
supply  options  the  company  is  considering  include: 

New  EPNG  Sales  Service.  EPNG  has  proposed  a  new 
sales  service  that  includes  a  monthly  demand  charge, 
a  commodity  charge  based  on  the  quantity  of  gas 
purchased,  and  a  gas  inventory  charge  (GIC).  A  GIC 
penalty  could  be  assessed  whenever  gas  purchases  fall 
below  a  target  level  tied  to  the  buyer’s  maximum  daily 
purchase  entitlement.  Southwest  Gas  can  choose  a 
new  daily  purchase  entitlement  between  zero  and  its 


454  /  Avery  et  al. 


RENO  JO 
INCLINE  VILLAGE  /  # 

carson  cnyjff$ 

STATELINE 
SOUTH  LAKE  TAHOE 


SOUTHWEST  GAS  CORPORATION 
EL  PASO  NATURAL  GAS 
PACIFIC  GAS  &  ELECTRIC 
NORTHWEST  PIPELINE 
PROPOSED  STORAGE  FACILITY 
COMMUNITIES  SERVED  BY  SOUTHWEST  GAS 
COMMUNITIES  SERVED  BY  OTHER  COMPANIES 
DIVISION  HEADQUARTERS 
LNG  FACILITY 
LPG  FACILITY 


Figure  1.  Southwest  Gas  Corporation’s  major  facilities  and  transmission  lines.  Also  shown  are  the  principal 
communities  where  the  corporation  supplies  gas  either  as  a  wholesaler  or  distributor. 
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current  entitlement  level.  FERC  action  is  required  on 
EPNG’s  tariff  proposal  before  the  exact  provisions  of 
this  service  will  be  known. 

Long-Term  Direct  Purchases.  Southwest  Gas  has 
received  proposals  from  gas  producers  for  long-term 
gas  sales  over  periods  of  up  to  five  years.  Contract 
provisions  typically  include  a  commodity  rate  tied  to 
an  index  of  spot  market  prices,  and  either  a  monthly 
demand  charge  or  a  minimum  purchase  requirement 
with  some  form  of  deficiency  penalty.  Gas  purchased 
directly  from  producers  would  be  transported  under 
firm  transportation  arrangements  with  EPNG. 

Pataya  Storage.  Southwest  Gas  is  considering  devel¬ 
opment  of  a  new  gas  storage  field  in  northwest 
Arizona.  Storage  would  allow  Southwest  Gas  to  reduce 
its  peak-day  entitlements  with  suppliers  (and  associ¬ 
ated  demand  charges)  and  purchase  additional  spot 
gas  in  the  off-peak  season.  The  storage  injection, 
inventory  capacity  and  withdrawal  sizing  have  been 
optimized  with  the  model  under  various  postulated 
future  market  conditions.  The  results  indicate  that 
building  the  storage  could  save  the  company’s  cus¬ 
tomers  $5  to  $40  million  per  year.  The  company  is 
closely  monitoring  market  conditions  which  affect  this 
analysis  and  when  forecasts  firm  up  the  company  will 
decide  whether  or  not  to  go  ahead  with  construction. 
If  the  decision  is  to  build,  the  model  will  be  used  to 
determine  optimum  injection,  size  and  withdrawal 
specifications.  The  Pataya  storage  could  become  avail¬ 
able  as  early  as  1993. 

The  Southwest  Gas  Supply  Planning  Group  is  cur¬ 
rently  testing  both  near-term  and  long-term  supply 
options  for  the  Southern  System.  One  analyst  is 
employed  full-time  to  make  model  runs  and  interpret 
results.  To  date,  over  50  gas  contract  options  and 
combinations  of  options  have  been  tested  over  one, 
three,  and  five-year  timeframes.  Promising  supply 
alternatives  are  then  tested  against  a  range  of  forecasts 
for  gas  sales  requirements  and  spot  price  conditions. 
A  typical  five-year  scenario  involves  4,600  constraints, 
12,900  variables  and  61,000  nonzero  coefficients. 
Optimization  of  each  scenario  takes  one  to  ten  min¬ 
utes  of  CPU  time  to  run  on  an  IBM  3090  operating 
under  MVS/XA. 


9.  A  CASE  STUDY  AT  QUESTAR 
CORPORATION 

Questar  Corporation  is  a  regional,  integrated  energy 
company  headquartered  in  Salt  Lake  City,  Utah 


Figure  2.  Questar  oil  and  gas  production  areas.  Much 
of  the  natural  gas  is  transported  by  Questar 
Pipeline  Company  to  Mountain  Fuel 
Supply  Company. 


(Questar  1989).  Questar  engages,  via  subsidiaries 
Wexpro/Celsius,  in  oil  and  gas  exploration  and  pro¬ 
duction  throughout  the  Rocky  Mountain  region  and 
elsewhere  (see  Figure  2).  Questar  Pipeline  operates  an 
interstate  pipeline  transmission  system  in  Colorado, 
Wyoming  and  Utah,  and  owns  underground  storage 
facilities  in  Wyoming  and  Utah  (see  Figure  3).  Questar 
Corporation  also  has  a  marketing  affiliate,  Questar 
Energy  Company,  which  pursues  unregulated  gas 
acquisition,  gathering,  compression  and  sales.  In  addi¬ 
tion,  Questar  operates  an  LDC,  Mountain  Fuel  Supply 
Company,  which  provides  retail  gas  distribution  to 
about  one-half  million  customers  in  Utah,  southwest¬ 
ern  Wyoming  and  southern  Idaho.  Thus,  Questar 
engages  in  virtually  all  aspects  of  the  natural  gas 
industry. 

Questar’s  three  principal  subsidiaries — engaged  in 
exploration  and  production,  storage,  transmission  and 
local  distribution — have  all  been  affected  profoundly 
by  deregulation.  Wexpro  wants  to  develop  and  pro¬ 
duce  more  gas  than  it  now  contributes  to  Questar 
customers.  Questar  Pipeline  seeks  to  provide  reliable 
and  profitable  transmission  and  storage  services,  and 
Mountain  Fuel,  as  an  LDC,  seeks  to  maximize  return 
on  stockholders’  equity  subject  to  existing  contractual 
commitments  with  suppliers  and  subject  to  approval 
from  regulatory  agencies.  These  goals  sometimes 
conflict. 

A  case  study  serves  to  illustrate  some  of  the  issues 
raised  by  deregulation.  Major  industrial  customers 
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Figure  3.  Questar  Pipeline  transmission  system.  Producers  and  other  pipelines  supply  natural  gas  which  is  stored 
or  distributed  to  the  local  distribution  company,  Mountain  Fuel. 


have  traditionally  contracted  for  inexpensive,  inter¬ 
ruptible  gas  from  Mountain  Fuel — interruptible  cus¬ 
tomers  provide  a  handy  way  to  balance  demand. 
However,  under  open  access,  most  of  these  customers 
have  been  able  to  contract  directly  with  gas  suppliers 
other  than  Mountain  Fuel,  and  have  obtained  inter¬ 
ruptible  transportation  from  Questar  Pipeline  and 
Mountain  Fuel.  If  the  other  gas  suppliers  have  trouble 
meeting  their  commitments  to  the  industrial  cus¬ 
tomers,  Mountain  Fuel  has  firm  and  interruptible 
concurrent  sales  agreements  with  many  customers  to 
make  up  the  deficit. 

The  evolution  from  a  single  source  of  gas  for  end- 
use  customers  to  a  plethora  of  providers  coupled  with 
accessible  transportation  has  had  adverse  conse¬ 
quences  for  Questar  Pipeline.  Questar  Pipeline  must 
negotiate  a  portfolio  of  supply  contracts  to  ensure 
long-term  availability  of  baseload  gas,  and  must  pro¬ 
vide  storage  and  peak-day  capacity  to  Mountain  Fuel 
for  gas  from  many  sources  of  varying  reliability.  In 
contrast,  Mountain  Fuel’s  major  industrial  customers 
who  have  converted  to  direct  producer  contracts  with 
firm  concurrent  fallback  have  essentially  become  part 
of  the  peak-load  problem  rather  than  contributing  to 


its  solution.  In  the  past  few  years,  the  proportion  of 
interruptible  demand  as  a  fraction  of  total  demand 
has  decreased  significantly. 

The  three  Questar  subsidiaries  decided  to  rationalize 
their  individual  competitive  strategies  in  the  context 
of  the  entire  gas  supply  system.  For  this  planning 
exercise,  a  nine-year  horizon  was  chosen  at  a  monthly 
level  of  detail.  Analysis  and  survey  of  industry  litera¬ 
ture  resulted  in  three  alternate  future  cost  and  demand 
load  factor  (the  proportion  of  maximum  contracted 
gas  which  is  take-or-pay)  scenarios:  nominal  price  and 
demand  increases,  moderate  increases  in  price  and 
demand,  and  a  scenario  with  a  moderate  increase  in 
demand  and  a  sharp  increase  after  a  few  years.  (The 
last  scenario,  in  which  a  sharp  price  increase  occurs, 
is  a  “gas  bubble”  scenario  which  many  industry  ana¬ 
lysts  predict.  Evaluation  of  such  a  scenario  makes 
good  business  sense.) 

With  these  scenarios,  and  the  existing  portfolio  of 
gas  contracts,  Questar  Pipeline  provided  modeling 
support  to  render  a  viable  joint  policy.  The  key  goals 
were  providing  projected  reliable  gas  supplies  at 
reasonable  prices  while  balancing  the  requirements  of 
both  customers  and  subsidiaries.  Fifty  to  sixty  model 
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excursions  were  evaluated,  each  with  about  6,200 
constraints,  12,000  variables  and  8 1 ,000  nonzero  coef¬ 
ficients.  Each  scenario  ran  in  9  megabytes  (2  mega¬ 
bytes  for  the  decision  support  shell,  2  megabytes  for 
optimization,  and  the  remainder  for  operating  system 
overhead)  for  about  3  minutes  on  an  IBM  3084-Q 
running  under  MVS/XA-SP. 

Surprisingly,  the  solutions  indicated  that  a  modest 
amount  of  additional  storage  and  peak-day  capacity 
will  provide  a  great  deal  of  leverage  to  defer  and  reduce 
the  magnitude  of  contract  purchases  for  the  purpose 
of  meeting  future  peak  demands.  In  particular,  the 
solutions  involve  expanded  use  of  the  Clay  Basin 
storage  field  (see  Figure  3),  which  is  a  depleted  “cap- 
rock”  gas  field  in  northeastern  Utah  with  7.5  BCF 
base  and  125  MMCF  peak  capacity  dedicated  to 
Questar  Pipeline.  This  can  be  augmented  at  attractive 
cost  by  a  2-3  BCF  base  and  50  MMCF  peak  capacity, 
thus  deferring  and  reducing  the  need  for  additional 
purchase  agreements  over  the  next  several  years. 
While  this  conclusion  may  not  appear  controversial, 
many  competing  proposals  of  varying  cost  and  com¬ 
plexity  required  evaluation  to  the  satisfaction  of  com¬ 
peting  executive  constituents.  It  is  in  this  environment 
that  good  models,  especially  good  optimization 
models,  excel. 

10.  CONCLUSIONS 

This  model  formulation  has  several  important  advan¬ 
tages  over  those  restricted  to  generalized  networks 
(e.g.,  McBride  1986,  Planmetrics  1988).  First,  a  com¬ 
plete  representation  of  the  supply  planning  problem 
is  provided.  At  the  strategic  level,  the  relationships 
among  gas  purchase  contracts,  storage,  and  transpor¬ 
tation  in  a  gas  supply  portfolio  are  captured.  At  the 
dispatch  level,  the  model  includes  important  operating 
and  contractual  constraints.  These  include  the  system 
network,  contract  pricing  terms,  and  the  relationships 
between  inventory  levels  and  injection  and  withdrawal 
capabilities.  The  diversity  of  available  contractual 
terms  mandates  incorporation  of  many  complicating 
constraints  and  variables.  That  the  model  accurately 
incorporates  the  right  amount  of  detail  is  demon¬ 
strated  by  the  model’s  wide  acceptance  in  the  natural 
gas  industry  (Rosenkranz  1989). 

From  the  perspective  of  the  gas  supply  planner,  it  is 
also  important  that  the  problem  is  actually  solved 
optimally.  In  contrast,  it  has  been  fashionable  to  sug¬ 
gest  Lagrangian  relaxation  techniques  to  render  easier 
problems  (i.e.,  generalized  networks,  or  even  pure 
networks  assuming  no  transmission  losses)  and 
attempt  to  satisfy  remaining  complicating  constraints 


heuristically,  i.e.,  with  no  assurance  of  success.  The 
Lagrangian  approach  allows  the  use  of  much  simpler 
optimization  packages,  e.g.,  Bradley,  Brown  and 
Graves  (1977),  and  Brown  and  McBride  (1984).  How¬ 
ever,  the  complicating  constraints  in  natural  gas  con¬ 
tracts  constitute  the  core  of  the  business  aspects  of  the 
problem — indirect  treatment  of  these  constraints 
inevitably  admits  infeasible  and  suboptimal  solutions 
to  the  problem  at  hand.  At  best,  this  greatly  compli¬ 
cates  comparisons  of  competing  scenarios — at  worst, 
comparisons  are  meaningless.  The  gas  supply  planner 
is  well  advised  to  depend  upon  complete,  feasible, 
optimal  analyses  of  complex  scenarios  for  presenta¬ 
tion  before  company  management  and  regulatory 
agencies. 

NOTES 

1  The  U.S.  gas  industry  commonly  uses  BTUs, 
THERMs  (10s  BTUs),  MMBTUs  (10s  BTUs),  cubic 
feet  (denoted  cf,  approximately  103  BTU),  MCF 
(103  cf),  MMCF  (106  cf),  and  BCF  (109  cf).  A  BTU 
may  better  be  expressed  outside  the  U.S.  as  1,055 
joules  or  252  gram  calories. 

APPENDIX 

Subscripts 

k  nodes; 

Jj  the  set  of  child  nodes  associated  with  parent 
node;; 

jk  the  transportation  agreement  from  node  j  to 
node  k ; 

t,  t  time  period; 

F  the  last  time  period  in  the  planning  horizon; 

/  side  constraint; 

S  season,  i.e.,  a  contiguous  set  of  periods; 

S  the  first  time  period  of  season  S'; 

T  a  contract  period,  a  contiguous  set  of  periods; 

T  the  first  time  period  of  contract  period  T ; 

T  the  last  time  period  of  contract  period  T ; 

Tt  the  set  of  periods  associated  with  side 
constraint  /; 

L  the  set  of  side  constraints; 

A  the  set  of  transportation  agreement  pairs  jk. 

Constants 

DAYS,  the  number  of  days  in  time  period  t\ 
DMDMj,  the  demand  at  node  j  during  time 
period  t\ 

LOADj,  peak-day  load  factor  at  demand  node  j 
during  time  period  t; 
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MNDWj,  minimum  daily  withdrawal  from  stor¬ 
age  node  j  during  time  period  t; 

MXDWj,  maximum  daily  withdrawal  from  stor¬ 
age  node  j  during  time  period  t; 

MNDIji  minimum  daily  injection  at  storage 
node  j  during  time  period  t; 

MXDIjt  maximum  daily  injection  at  storage 
node  j  during  time  period  t; 

TINVj,  minimum  (target)  inventory  level  at 
storage  node  j  at  end  of  time  period  t; 

SCAPji  maximum  inventory  level  at  storage 
node  j  at  end  of  period  f, 

MNCP/t  minimum  allowable  capacity  at  poten¬ 
tial  storage  node  j  over  contract  period 
T; 

MXCPjt  maximum  allowable  capacity  at  poten¬ 
tial  storage  node  j  over  contract  period 
T; 

MNMWjt  minimum  maximum  daily  withdrawal 
from  potential  storage  node  j  over  con¬ 
tract  period  T ; 

MXMWjt  maximum  maximum  daily  withdrawal 
from  potential  storage  node  j  over  con¬ 
tract  period  T ; 

MXDPj ,  maximum  daily  purchase  from  pipe¬ 
line  contract  node  j  during  time  period 

t; 

MNDPji  minimum  daily  purchase  from  pipeline 
contract  node  j  during  time  period  t; 

MXSPjs  maximum  seasonal  purchase  from 

pipeline  contract  node  j  during  season 

MDPFjs  minimum  daily  purchase  fraction  from 
pipeline  contract  node  j  over  contract 
period  T; 

MNPFjs  minimum  purchase  fraction  from  pipe¬ 
line  contract  node  j  in  season  S; 

MNPPjs  minimum  purchase  penalty  for  pipe¬ 
line  contract  node  j  in  season  S; 

MXMDjt  maximum  allowable  maximum  daily 
purchase  from  potential  pipeline  con¬ 
tract  node  j  over  contract  period  T ; 

SDA  Yjt  the  number  of  days  in  a  season  for 
potential  contract  period  T  if  mini¬ 
mum  seasonal  purchase  is  based  on 
maximum  daily  purchase  ujT  (con¬ 
straint  ( MNPF ')  is  used);  otherwise, 
SDA  YjT  =£  0  indicates  that  constraint 
(MNPF)  is  used,  i.e.,  minimum  sea¬ 
sonal  purchase  is  based  on  maximum 
seasonal  purchase  zJT; 

JCAPjk,  the  joint  firm  and  interruptible  daily 


capacity  on  potential  arc  jk  during 
period  t ; 

MXMFjk  maximum  firm  daily  capacity  on 
potential  arc  jk-, 

MXMIjk  maximum  interruptible  daily  capacity 
on  potential  arc  jk ; 

MNDTjk,  minimum  daily  throughput  on  existing 
arc  jk  during  time  period  t; 

MXDTjkt  maximum  daily  throughput  on  existing 
arc  jk  during  time  period  t; 

RELFj,  reliability  factor  for  node  j  during  time 
period  t; 

RELFjk,  reliability  factor  for  firm  arc  jk  during 
time  period  t; 

RELIjk,  reliability  factor  for  interruptible  arc  jk 
during  time  period  t; 

INTAjki  interruptibility  factor  for  interruptible 
arc  jk  during  time  period  t; 

JMXMjt  maximum  joint  maximum  daily  pur¬ 
chase  on  potential  joint  contract  node  j 
over  contract  period  T ; 

FCMXjt  maximum  fraction  of  joint  maximum 
daily  purchase  children  of  parent  j  may 
have  over  contract  period  T ; 

LHSV/  the  left-hand  side  value  for  side 
constraint  /; 

RHSVi  the  right-hand  side  value  for  side 
constraint  /; 

mjk  the  flow  multiplier  (1  -  loss)  on  arc  jk; 
ijo  the  initial  inventory  at  storage  node  j; 
mj  the  flow  multiplier  on  injection  into 
storage  node  j ; 

mj  the  flow  multiplier  on  withdrawal  from 
storage  node  j; 

djT  the  lower  maximum  daily  purchase/ 
maximum  seasonal  purchase  constant 
for  node  j  over  contract  period  T ; 
dyr  the  upper  maximum  daily  purchase/ 
maximum  seasonal  purchase  constant 
for  node  j  over  contract  period  T ; 
djT  the  joint  contract  lower  maximum 
daily  purchase/maximum  seasonal 
purchase  constant  seen  at  parent  node  j 
over  contract  period  T ; 
djT  the  joint  contract  upper  maximum 
daily  purchase/maximum  seasonal 
purchase  constant  seen  at  parent  node  j 
over  contract  period  T ; 
fj  the  daily  fraction  of  inventory  at  stor¬ 
age  node  j  which  can  be  withdrawn 
(withdrawal  function  slope); 
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hj  the  least  daily  withdrawal  amount  at 
storage  node  j  (withdrawal  function  y- 
axis  intercept); 

f,  the  daily  fraction  of  inventory  at  stor¬ 
age  node  j  which  can  be  injected  (nega¬ 
tive  of  injection  function  slope); 
hj  the  maximum  daily  injection  amount 
at  storage  node  j  given  that  MXDI  is 
not  binding  (injection  function  y-axis 
intercept)  at  storage  node  j; 
r,  the  discount  factor  for  period  V, 
ajk/  the  coefficient  associated  with  arc  jk  in 
side  constraint  /; 

CCOMj,  the  cost  of  commodity  at  supply  node  j 
during  time  period  t; 

CTFMjk,  the  cost  of  firm  transportation  on 

potential  arc  jk  during  time  period  f, 

CTINjk,  the  cost  of  interruptible  transportation 
on  potential  arc  jk  during  time  period 
t; 

CDFMjkr  the  demand  charge  for  firm  transporta¬ 
tion  on  potential  arc  jk,  contract  period 
T; 

CDINjt  the  demand  charge  for  interruptible 

transportation  on  potential  arc  jk  over 
contract  period  T ; 

CTRNjk,  the  cost  of  transportation  on  existing 
arc  jk  during  time  period  t\ 

CINJj  the  cost  of  injections  at  storage  node  j; 

CINVj  the  daily  inventory  holding  cost  at  stor¬ 
age  node  j; 

CWDRj  the  withdrawal  charge  at  storage  node 

j\ 

CCAPjt  the  cost  of  capacity  at  potential  storage 
node  j  over  contract  period  T\ 

CMDWjt  the  cost  of  maximum  daily  withdrawal 
at  potential  storage  node  j  over  con¬ 
tract  period  T ; 

CDMDjt  the  demand  charge  at  potential  supply 
node  j  over  period  T ; 

CDMSjt  the  seasonal  demand  charge  at  poten¬ 
tial  supply  node  j  over  period  T ; 

CJDDjr  the  demand  charge  for  joint  contract 
seen  at  parent  potential  supply  node  j 
over  contract  period  T ; 

CJDSjt  the  seasonal  demand  charge  for  joint 

contract  seen  at  parent  potential  supply 
node  j  for  period  T ; 

CPENj  the  penalty  charge  for  not  meeting 
demand  at  node  j; 

CLPN/  the  penalty  for  violating  left-hand  side 
of  side  constraint  /; 


CRPN/  the  penalty  for  violating  right-hand  side 
of  side  constraint  /. 

Variables 

x'Jkt  the  flow  from  node  j  to  node  k  in  time  period  t 
on  potential  firm  arc; 

xjk,  the  flow  from  node  j  to  node  k  in  time  period  t 
on  potential  interruptible  arc; 

Xjk,  the  flow  on  existing  arc  jk  in  time  period  t; 
xjk,  the  total  flow  on  all  arcs  from  node  j  to  node  k 
in  time  period  t; 

xf,  the  peak-day  flow  from  node  j  to  node  k  in  time 
period  t  on  potential  firm  arc; 
xjk,  the  peak-day  flow  from  node  j  to  node  k  in  time 
period  t  on  potential  interruptible  arc; 

Xjk,  the  peak-day  flow  on  existing  arc  jk  in  time 
period  t; 

Xjk,  the  total  peak-day  flow  on  all  arcs  from  node  j 
to  node  k  in  time  period  t; 
vj,  injections  at  storage  node  j  in  time  period  t; 

Wj,  withdrawals  from  storage  node  j  in  time 
period  P, 

fr  the  variable  capacity  on  potential  storage  node  j, 
contract  period  T; 

WjT  the  variable  upper  bound  on  withdrawals  from 
node  j;  contract  period  T; 

bj,  the  flow  directly  bypassing  storage  node  j  in  time 
period  t; 

Wj,  withdrawals  from  storage  node  j  on  peak  day  of 
time  period  t\ 

yjkT  the  variable  upper  bound  on  firm  capacity  in  arc 
jk,  contract  period  T ; 

yjkr  the  variable  upper  bound  on  interruptible  trans¬ 
portation  capacity  in  arc  jk,  contract  period  T ; 
ujT  the  variable  maximum  daily  purchase  for  poten¬ 
tial  pipeline  contract  node  j,  contract  period  T; 
zjT  the  variable  maximum  seasonal  purchase  for 
potential  pipeline  contract  node  j,  contract 
period  T; 

Ujr  the  variable  maximum  daily  purchase  for  poten¬ 
tial  joint  pipeline  contract  whose  parent  is  node 
j,  contract  period  T ; 

ZjT  the  variable  maximum  seasonal  purchase  for 
potential  joint  pipeline  contract  whose  parent  is 
node  j,  contract  period  T ; 
sf,  the  unmet  demand  at  node  j  in  time  period  t; 
sjj  the  violation  of  minimum  purchase  (minimum 
bill)  for  contract  node  j  for  season  S; 
sj  the  unmet  peak-day  demand  at  node  j  in  time 
period  t. 
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Period-by-Period  Constraints 

Spot  purchase  node  j: 

MNDPj.DAYS ,  =£  Z  xJkt  =£  MXDP]tDAYS, 

k 

for  all  t.  (. MXMP ) 

Existing  pipeline  contract  node  j: 

MNDPjtDA  YSt  =£  Z  xJk,  «  MXDPjtDA  YS, 

k 

for  all  /  t MXMP ) 

-s%  +  MXSPjsMNPFjs  Z  I  x,fe  =£  MXSPjs 

k  ,GS 

for  all  S.  ( MXSP ) 

Potential  pipeline  contract  node  j: 

Z  Xyfa  -  DA  YS,  u,T  «  0 

k 

for  all  t  G  T,  for  all  T  ( MXMM ) 
0  «  Z  xjkt  -  MDPFj,  DA  YSt  uJT 

k 

for  all  t  G  T,  for  all  T  (. MNMP ) 
Xk>  -  MNPFjsZjt  for  all  S  G  T, 

k  teS 

for  all  T  if  SDA  Yjr  «  0  (MNPF) 
-sM  «  Z  Z  x>fa  -  SDAYjtMNPFjsUjt 

k  t£S 

for  all  5  G  T,  for  all  T  if  SDA  YJ  f  >  0  {MNPF’) 
Z  Z  x#r  -  ZjT  o 

k  keS 

for  all  SET,  for  all  T  (SEAS) 
0  =£  zjT  -  djTUjT  for  all  T  (DSLO) 

ZjT  -  HJTUjT  «  0  for  all  T  (DSUP) 

MNDPjr  *=  ujT  =5  MXMDjt,  0  ZjT  for  all  T. 

Z  Z  Kj’kt  +  Z  Xjkt  -  DAYS,  UjT  =S  0 

j'EiJj  k  k 

for  all  t  G  T,  for  all  T  (JMXM) 

Z  Z  Z  x/'fa  +  Z  Z  *jkt  -  zjt  --=  o 

j'eJj  k  te.S  k  teS 

for  all  .S’  G  T  for  all  T  (JSEA) 

Z  Z  xj’kt  ~  DA  YS/FCMXjtUjt  0 

j’ejj  k 

for  all  t  G  T,  for  all  T  (JCPM) 


0  s=  zjt  ~  djTUjr  for  all  T  (JDSL) 

ZjT  -  cljiijr  0  for  all  T  (. JDSU ) 

0  UjT  =£  JMXMjt,  0  =£  z]T  for  all  T. 

Existing  storage  node  j: 

1  h 

fjDA  YS,  Wjtr  +  ^  ~m'Vjr  +  ^  Wjr  *  ij0  +  ij 

for  all  t  ( WINY ) 

1  h' 

Jj  DAYS \  Vj‘  +  ?,  m'VjT  +  ?,  ~Wjr  ^  ~ij0  +  jf 

for  all  t  ( IJIV ) 

X  -mkjxkj,  +  bj,  +  Vj,  =  0  for  all  t  ( INPT ) 

k 

t  t 

-ijo  +  TINVj,  =s  Z  mjVjr  +  Z  ~wjr  SCAPj  -  ij0 

T  —  1  T=1 

for  all/  (STOR) 
-mJWj,  -  bjt  +  Z  Xjkt  =  0  for  all  t.  ( OUTP ) 

k 

MNDIjiDAYS,  Vj,  =s  MXDIjtDA YS„ 

DA  YS,MND Wjt  =s  wJt  =£  MXDWJtDAYS,  for  all  t. 
Potential  storage  node  j: 

W1NV,  IJIV,  INPT,  STOR,  OUTP,  and 

-ijo  +  TINVj,  ^  Z  mvjVjr  +  Z  ~Wjr 

T=  1  T  —  1 

for  all  t  (INVT) 

Wjt  DA  YSt  WjT  ^  0 

for  all  t  G  T,  for  all  T  ( VMDW ) 

MNCPjt  «  Tjr  «  MXCPjt  for  all  T 
MNMWjt  =£  wjT  MXMWjt  for  all  T 
0  <  Wj,  for  all  t 

MNDIjTDAYSt  =s  Vj,  =s  MXDIjrDAYS, 

for  all  t  G  T,  for  all  T. 

Demand  node  j: 

z  -rrikjKkj,  +  z  Xjkt  -DMDMj,  +  sf, 

k  k 

for  all  t.  (DMDM) 
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Junction  node  j: 

2  -mkjxkj,  +  2  xJkl  =  0  for  all  t.  ( JUNC ) 

k  k 

Existing  transportation  agreement  jk: 

MNDTjk,DA  YS,  <  xjk, MXDTjk,DA  YS ,  for  all  t. 
Potential  transportation  agreement  jk: 
xjkl  +  x;kt  JCAPjkrDA  YS, 

for  all  t  G  T,  for  all  T  ( JCAP ) 

x’kl  -  DAYS.yfr  =S  0 

for  all  t  G  T,  for  all  T  ( MXMF ) 

x"lkl  -  DAYS,INTAjklyfr  *£  0 

for  all  t  G  T,  for  all  T  ( MXMI ) 

0  =£  y jkT  MXMFjkT,  0  *£  yjkT  «  MXMIjkT  for  all  T 
0  «  xjkt  *£  DAYS, MXMFjkT,  0  «  xjkt 

DAYS,MXMIjkTINTAjkt  for  all  t  G  T,  for  all  T. 
Side  constraints  on  fixed  arcs: 

-sj  +  LHSV,  «  ^  2  ajk/Xjkl 

1 67)  U.k)eA 

RHSVi  +  st  for  all  /  G  L,  for  all  L.  ( ALNK ) 

Peak-Day  Constraints 

Demand  node  j: 

2  mkjxkJt  "f"  2  Nkt 

k  k 

-DMDMj,LOADj,/DA  YS,  +  sf, 

for  all  t.  ( PDMD ) 

Simple  junction  node  j: 

2  mkjxkj,  +  2  xjk,  =  0  for  all  t.  ( PJNC ) 

k  k 

Spot  purchase  node  j  or  existing  pipeline  contract 
node  j: 

0  D  xJk,  =£  MXDPj,  RELFj,  for  all  t.  ( PMDP ) 

k 

Potential  pipeline  contract  node  j: 

2  kjk,  -  RELFj,  ujT  =£  0 

k 

for  all  t  G  T,  for  all  T  ( PMMD ) 


2  2  Xj'fa  +  2  xjk,  -  RELFjtUjT  =£  0 

k  k 

for  all  t  G  F,  for  all  T  {PJMM) 

2  2  N’k,  -  FCMXjtUjt  0 

j’ejj  k 

for  all  t  G  T,  for  all  T.  ( PJCC ) 
Existing  storage  node  j: 

Yu-mkj  xkj,  -  mj  Wj,  +  2  ^kt  =  0  for  all  t  ( PSTR ) 

k  k 

\  1  1  fo. 

7  Wj,  +  2  -  mj  vjT  +  2  *f/T  =S  ij0  +  -jr 

Jjt  T=l  T=  1  Jjt 

for  all  l  ( P1NV) 

t  t 

wj,  +  2  ~mj  vjr  +  2  wjr  iJ0  for  all  t  (PTIV) 

T= 1  T= 1 

0  Wj,  *£  RELFj,  MXDWj,  for  all  t. 

Potential  storage  node  j: 

PSTR,  PINV,  PTIV,  and 

Wj,  -  RELFJTWjT  =s  0 

for  all  t  G  T,  for  all  T  (PVMW) 

0  =s  Wj,  for  all  t. 

Existing  transportation  agreement: 

MNDTjk,  xjk,  MXDTjk,RELFjk,  for  all  t. 
Potential  transportation  agreement  jk: 
xjk,  +  Xjk,  «  7C4P*  for  all  t  ( PTCP ) 

xf,  -  RELFjktyjkT  =s  0 

for  all  t  G  T,  for  all  7'  (PMXLj 

Xjk,  ~  RELIjk,y"kr  ^  0 

for  all  t  G  T,  for  all  T  ( PMXI ) 

0  =s  x,'c(,  0  s;  x'jk,  for  all  t. 

Costs 

xjk,  r,CTFMjk„  or  r,CTFMjk,  +  r,CCOMjt  if  j  is  a 
supply  node; 

x"k,  r,CTINjkt,  or  r,CTINJk,  +  r,CCOMjt  if  j  is  a 
supply  node; 
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Xjk,  r,CTRNjk,,  or  r,CTRNjk,  +  r,CCOMJt  if  j  is  a 
supply  node; 

Xjk,  r,CTFMjkl,  or  rtCTFMjk,  +  r,CCOMjt  if  j  is  a 
supply  node; 

x":,  r,CTINjk„  or  r,CTINjk,  +  r,CCOMjt  if  j  is  a 
supply  node; 

xik,  r,CTRNjk„  or  r,CTRNJk,  +  rtCCOMjt  if  j  is  a 
supply  node; 

vjt  r, CINJj  +  ZL,  m'j  r7 DA  YS\  CINV,; 
wj,  r, CWDRj  -  ZL,  rTDAYSTCINV }; 

ijr  rk  CCAPj/j 
WjT  yjjCMDWjT'i 
hj,  no  cost; 
wjt  r.CWDRj- 

yjkr  ZI=t  y,CDFMjkr\ 

y'JkT  ZlrnCDlNjkT- 
Ujr  E  r  ri  C DMDjj ; 

ZjT  ZhnCDMSjr ; 

Ujr  ztr  r,CJDDjT', 

Zjt  Zti  nCJDSjr, 
sf,  r,CPENj\ 
sfs  rs  CDMDjMNPPjs, 
sf,  r,CPENj- 

sj  CLPN/  times  appropriate  discount  rate; 
sf  CRPNi  times  appropriate  discount  rate. 
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